Correlation functions for time-dependent calculation of linear-response functions 
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We emphasize the importance of choosing an appropriate 
correlation function to reduce numerical errors in calculating 
the linear-response function as a Fourier transformation of a 
time-dependent correlation function. As an example we take 
dielectric functions of silicon crystal calculated with a time- 
dependent method proposed by Iitaka et al. [Phys. Rev. E 
56, 1222 (1997)]. 

PACS: 02.70.-c, 71.10.Fd, 85.30.Vw, 72.15.-v 

Recently, Iitaka it et al. proposed a linear scaling algo- 
rithm for calculating linear- response functions of nonin- 
teracting electrons . This algorithm has been success- 
fully applied to the calculation of dielectric functions of 
Si crystal 0] and Si nano crystallites They calculate 
the linear-response of B due to the perturbation Ae~ %ult 



X ba(u + if,) = ^ ^ dt e+^+ l ^K(t)jj , (1) 

where the double angular brackets indicate the statistical 
average over random vectors |<J>} and K (t) is the time- 
dependent correlation function defined by 

K{t) = 2 lm{<5>\6(E } - H)e +lHt Be- lHt 6{H - E f )A\$). 

(2) 

In the case of the dielectric function ep. a = 1 + knxfia 
with a, = x,y, z the perturbation and the response are 
the external electric field and polarization, respectively, 



(j\A\i) = (j\x a \i), 
(W) = ^(i\x0\j), 



(3) 
(4) 



where V = L 3 is the volume of the supercell. 

In Ref. H Iitaka et al. further modified the dipole 
moment operator x to the current operator j by partial 
integration 




dt e-" 4 , 

+ IT,) 



K 3 (t) 



V 



■Im(<Z>\0(E f -H)e 



+iHt 



-iHt 



0(E cut -H)0(H-E f )j a \$). 
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and calculated the dielectric function of silicon crystal 
consisting of 2 15 Si atoms in a cubic supercell of 16 3 
unit cells with energy resolution r, = 0.05 (eV). Each 
unit cell is divided into 8 3 cubic meshes. Iitaka et al. Q| 
used the empirical local pseudopotcntial developed by by 
Wang and Zunger in Ref. |!|] and one random vector for 
the statistical average. Figure |l| shows the result cal- 
culated with the same condition as in Ref. §. Iitaka 
et al. attributed the divergence of calculated response 
functions at low frequencies to the artificially introduced 
finite imaginary part r, of the frequency. 

However, closer examination suggests that the di- 
vergence may originate from the modification of the 
position-position correlation function to the current- 
current correlation function. Assuming that the dielec- 
tric function is analytic at zero frequency, we easily no- 
tice that the low-frequency component of Kj(t) has a 
frequency dependence of u 2- * 15 ". Therefore, numerical 
results calculated with Kj(t) may easily suffer round- 
off errors when the small amplitude of the low-frequency 
components is extracted from the large background of 
higher- frequency components. 

To check this scenario, we calculate the dielectric func- 
tion without modification of dipole operators to cur- 
rent operators. Instead we modified Eqs. (||) and (jij) 
to a finite-wave-number form to conform with periodic 
boundary conditions 
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where q is the smallest wave number in the supercell 
q = 2n/L. The results are shown in Fig. ^, where other 
conditions are the same as in Fig. [j]. Here we obtained 
the dielectric function without any divergence at low fre- 
quencies, while the fluctuations at high frequencies are 
larger. We think that this result confirms our scenario 
that the divergence originates from the use of Kj(t) in- 
stead of K(t). 

In summary, we found the importance of choosing an 
appropriate form of the correlation functions in comput- 
ing linear-response functions as a Fourier transformation 
of time-dependent correlation functions. By choosing 
the current-current correlation function for high frequen- 
cies and the position-position correlation function for low 
frequencies, we may accurately calculate linear-response 
functions at any frequency. 
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FIG. 1. e xx (ui) of silicon crystal calculated with cur- 
rent-current correlation function Kj(t). The system consists 
of 2 15 Si atoms in a cubic supercell of 16 3 unit cells. Each unit 
cell is divided into 8 3 cubic meshes. The energy resolution is 
r\ — 0.05(eV). We used the empirical local pseudopotential in 
Ref. M. (a) Real part and (b) imaginary part. 



FIG. 2. Same as Fig. |lj but calculated with posi- 
tion-position correlation function K(t). 
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Fig. 1a T.litaka 




Fig. 1b T.litaka 




Fig. 2a T.litaka 




Fig. 2b T.litaka 



